Comparative transcriptome analysis provides molecular insights into heterosis of waterlogging tolerance in Chrysanthemum indicum.

BACKGROUND: Heterosis breeding is one of the most important breeding methods for chrysanthemum. To date, the genetic mechanisms of heterosis for waterlogging tolerance in chrysanthemum are still unclear. This study aims to analyze the expression profiles and potential heterosis-related genes of two hybrid lines and their parents with extreme differences in waterlogging tolerance under control and waterlogging stress conditions by RNA-seq. RESULTS: A population of 140 F1 progeny derived from Chrysanthemum indicum (Nanchang) (waterlogging-tolerant) and Chrysanthemum indicum (Nanjing) (waterlogging-sensitive) was used to characterize the extent of genetic variation in terms of seven waterlogging tolerance-related traits across two years. Lines 98 and 95, respectively displaying positive and negative overdominance heterosis for the waterlogging tolerance traits together with their parents under control and waterlogging stress conditions, were used for RNA-seq. In consequence, the maximal number of differentially expressed genes (DEGs) occurred in line 98. Gene ontology (GO) enrichment analysis revealed multiple stress-related biological processes for the common up-regulated genes. Line 98 had a significant increase in non-additive genes under waterlogging stress, with transgressive up-regulation and paternal-expression dominant patterns being the major gene expression profiles. Further, GO analysis identified 55 and 95 transgressive up-regulation genes that overlapped with the up-regulated genes shared by two parents in terms of responses to stress and stimulus, respectively. 6,640 genes in total displaying maternal-expression dominance patterns were observed in line 95. In addition, 16 key candidate genes, including SAP12, DOX1, and ERF017 which might be of significant importance for the formation of waterlogging tolerance heterosis in line 98, were highlighted. CONCLUSION: The current study provides a comprehensive overview of the root transcriptomes among F1 hybrids and their parents under waterlogging stress. These findings lay the foundation for further studies on molecular mechanisms underlying chrysanthemum heterosis on waterlogging tolerance.

Multi-locus genome-wide association studies reveal the dynamic genetic architecture of flowering time in chrysanthemum.

KEY MESSAGE: The dynamic genetic architecture of flowering time in chrysanthemum was elucidated by GWAS. Thirty-six known genes and 14 candidate genes were identified around the stable QTNs and QEIs, among which ERF-1 was highlighted. Flowering time (FT) adaptation is one of the major breeding goals in chrysanthemum, a multipurpose ornamental plant. In order to reveal the dynamic genetic architecture of FT in chrysanthemum, phenotype investigation of ten FT-related traits was conducted on 169 entries in 2 environments. The broad-sense heritability of five non-conditional FT traits, i.e., budding (FBD), visible coloring (VC), early opening (EO), full-bloom (OF) and decay period (DP), ranged from 56.93 to 84.26%, which were higher than that of the five derived conditional FT traits (38.51-75.13%). The phenotypic variation coefficients of OF_EO and DP_OF were relatively large ranging from 30.59 to 36.17%. Based on 375,865 SNPs, the compressed variance component mixed linear model 3VmrMLM was applied for a multi-locus genome-wide association study (GWAS). As a result, 313 quantitative trait nucleotides (QTNs) were identified for the non-conditional FT traits in single-environment analysis, while 119 QTNs and 67 QTN-by-environment interactions (QEIs) were identified in multi-environment analysis. As for the conditional traits, 343 QTNs were detected in single-environment analysis, and 119 QTNs and 83 QEIs were identified in multi- environment analysis. Among the genes around stable QTNs and QEIs, 36 were orthologs of known FT genes in Arabidopsis and other plants; 14 candidates were mined by combining the transcriptomics data and functional annotation, including ERF-1, ACA10, and FOP1. Furthermore, the haplotype analysis of ERF-1 revealed six elite accessions with extreme FBD. Our findings contribute to the understanding of dynamic genetic architecture of FT and provide valuable resources for future chrysanthemum molecular breeding programs.

Genetic architecture and genomic prediction of plant height-related traits in chrysanthemum.

Plant height (PH) is a crucial trait determining plant architecture in chrysanthemum. To better understand the genetic basis of PH, we investigated the variations of PH, internode number (IN), internode length (IL), and stem diameter (SD) in a panel of 200 cut chrysanthemum accessions. Based on 330 710 high-quality SNPs generated by genotyping by sequencing, a total of 42 associations were identified via a genome-wide association study (GWAS), and 16 genomic regions covering 2.57 Mb of the whole genome were detected through selective sweep analysis. In addition, two SNPs, Chr1_339370594 and Chr18_230810045, respectively associated with PH and SD, overlapped with the selective sweep regions from FST and π ratios. Moreover, candidate genes involved in hormones, growth, transcriptional regulation, and metabolic processes were highlighted based on the annotation of homologous genes in Arabidopsis and transcriptomes in chrysanthemum. Finally, genomic selection for four PH-related traits was performed using a ridge regression best linear unbiased predictor model (rrBLUP) and six marker sets. The marker set constituting the top 1000 most significant SNPs identified via GWAS showed higher predictabilities for the four PH-related traits, ranging from 0.94 to 0.97. These findings improve our knowledge of the genetic basis of PH and provide valuable markers that could be applied in chrysanthemum genomic selection breeding programs.

A group VIIIa ethylene-responsive factor, CmERF4, negatively regulates waterlogging tolerance in chrysanthemum.

Ethylene-responsive factors (ERF) play an important role in plant responses to waterlogging stress. However, the function and mechanism of action of ERFVIII in response to waterlogging stress remain poorly understood. In this study, we found that expression of the ERF VIIIa gene CmERF4 in chrysanthemum was induced by waterlogging stress. CmERF4 localized to the nucleus when expressed in tobacco leaves. Yeast two-hybrid and luciferase assays showed that CmERF4 is a transcriptional inhibitor. CmERF4 overexpression in chrysanthemum reduced plant waterlogging tolerance, whereas overexpression of the chimeric activator CmERF4-VP64 reversed its transcriptional activity, promoting higher waterlogging tolerance than that observed in wild-type plants, indicating that CmERF4 negatively regulates waterlogging tolerance. Transcriptome profiling showed that energy metabolism and reactive oxygen species (ROS) pathway-associated genes were differentially expressed between CmERF4-VP64 and wild-type plants. RT-qPCR analysis of selected energy metabolism and reactive oxygen species-related genes showed that the gene expression patterns were consistent with the expression levels obtained from RNA-seq analysis. Overall, we identified new functions of CmERF4 in negatively regulating chrysanthemum waterlogging tolerance by modulating energy metabolism and ROS pathway genes.

Multi-locus genome-wide association study and genomic prediction for flowering time in chrysanthemum.

MAIN CONCLUSION: Multi-locus GWAS detected several known and candidate genes responsible for flowering time in chrysanthemum. The associations could greatly increase the predictive ability of genome selection that accelerates the possible application of GS in chrysanthemum breeding. Timely flowering is critical for successful reproduction and determines the economic value for ornamental plants. To investigate the genetic architecture of flowering time in chrysanthemum, a multi-locus genome-wide association study (GWAS) was performed using a collection of 200 accessions and 330,710 single-nucleotide polymorphisms (SNPs) via 3VmrMLM method. Five flowering time traits including budding (FBD), visible colouring (VC), early opening (EO), full-bloom (OF) and senescing (SF) stages, plus five derived conditional traits were recorded in two environments. Extensive phenotypic variations were observed for these flowering time traits with coefficients of variation ranging from 6.42 to 38.27%, and their broad-sense heritability ranged from 71.47 to 96.78%. GWAS revealed 88 stable quantitative trait nucleotides (QTNs) and 93 QTN-by-environment interactions (QEIs) associated with flowering time traits, accounting for 0.50-8.01% and 0.30-10.42% of the phenotypic variation, respectively. Amongst the genes around these stable QTNs and QEIs, 21 and 10 were homologous to known flowering genes in Arabidopsis; 20 and 11 candidate genes were mined by combining the functional annotation and transcriptomics data, respectively, such as MYB55, FRIGIDA-like, WRKY75 and ANT. Furthermore, genomic selection (GS) was assessed using three models and seven unique marker datasets. We found the prediction accuracy (PA) using significant SNPs identified by GWAS under SVM model exhibited the best performance with PA ranging from 0.90 to 0.95. Our findings provide new insights into the dynamic genetic architecture of flowering time and the identified significant SNPs and candidate genes will accelerate the future molecular improvement of chrysanthemum.

FUL homologous gene CmFL1 is involved in regulating flowering time and floret numbers in Chrysanthemum morifolium.

Flowering time and floret numbers are important ornamental characteristics of chrysanthemums that control their adaptability and inflorescence morphology, respectively. The FRUITFULL (FUL) gene plays a key role in inducing flowering and inflorescence meristem development. In this study, we isolated a homolog of the MADS-box gene FUL, CmFUL-Like 1 (CmFL1), from chrysanthemum inflorescence buds. Quantitative RT-PCR and in situ analyses showed that CmFL1 was strongly expressed in young inflorescence buds. Overexpression of CmFL1 caused early flowering while co-suppression expression of CmFL1 increased the number of florets. Furthermore, the floral promoting factors CmSOC1, CmFDL1, and CmLFY were up-regulated in the shoot tips of transgenic plants. In addition, RNA-seq analysis of the transgenic plants suggested that certain differentially expressed genes (DEGs)-such as MADS-box, homeobox family, and ethylene pathway genes-may be involved in the inflorescence meristem development. GO pathway enrichment analysis showed that the differentially transcribed genes enriched the representation of the carbohydrate metabolic pathway, which is critical for flower development. Overall, our findings revealed the conserved function of CmFL1 in controlling flowering time along with a novel function in regulating the number of florets.

Analyses of a chromosome-scale genome assembly reveal the origin and evolution of cultivated chrysanthemum.

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a globally important ornamental plant with great economic, cultural, and symbolic value. However, research on chrysanthemum is challenging due to its complex genetic background. Here, we report a near-complete assembly and annotation for C. morifolium comprising 27 pseudochromosomes (8.15 Gb; scaffold N50 of 303.69 Mb). Comparative and evolutionary analyses reveal a whole-genome triplication (WGT) event shared by Chrysanthemum species approximately 6 million years ago (Mya) and the possible lineage-specific polyploidization of C. morifolium approximately 3 Mya. Multilevel evidence suggests that C. morifolium is likely a segmental allopolyploid. Furthermore, a combination of genomics and transcriptomics approaches demonstrate the C. morifolium genome can be used to identify genes underlying key ornamental traits. Phylogenetic analysis of CmCCD4a traces the flower colour breeding history of cultivated chrysanthemum. Genomic resources generated from this study could help to accelerate chrysanthemum genetic improvement.

CmERF5-CmRAP2.3 transcriptional cascade positively regulates waterlogging tolerance in Chrysanthemum morifolium.

Waterlogging stress affects plant growth by limiting root respiration and reducing yield and economic value. Therefore, identifying genes involved in regulating waterlogging stress is vital. This study reports the ethylene-responsive VII transcription factor (CmRAP2.3) in the chrysanthemum. Subcellular localization and transactivation assay analyses revealed that CmRAP2.3 was localized in the nucleus and possessed transactivation activity. Overexpression of CmRAP2.3 in chrysanthemum was found to enhance waterlogging tolerance by decreasing reactive oxygen species (ROS) levels. Furthermore, we found that the transcription factor CmERF5 binds to GCC-like motifs in the CmRAP2.3 promoter region and activates CmRAP2.3 expression. Additionally, CmERF5 overexpression maintained a low ROS level and improved chrysanthemum waterlogging tolerance. Taken together, this study shows a molecular mechanism by which CmERF5 transcriptionally activates CmRAP2.3 to reduce waterlogging stress via the ROS pathway in the chrysanthemum.

Evaluation of a Yeast Two-Hybrid Library by High-Throughput Sequencing.

Yeast two-hybridization (Y2H) is a classical method to study protein-protein interactions in organisms, and the top limitation of Y2H is the high ratio of false negatives and false positives. The most efficient way to reduce this error is to improve the quality of the library. The traditional library quality evaluation method can only inform us of the capacity of the library and a very limited number of library insert lengths. Therefore, we developed a new method to evaluate the library by using only a 10 ng library, amplifying the inserted fragments through 15 cycles of PCR, and then carrying out high-throughput sequencing. This method can eliminate the randomness and one-sidedness of the traditional method and can be used to obtain key indicators, such as the number of inserted genes and gene abundance, to effectively evaluate the quality of the library. In addition, the new library quality assessment method can also reveal the gene sequences of species. This method is expected to greatly accelerate PPI research on nonmodel species.

The core regulatory networks and hub genes regulating flower development in Chrysanthemum morifolium.

KEY MESSAGE: The study has facilitated important insights into the regulatory networks involved in flower development in chrysanthemum (Asteraceae), and is informative with respect to the mechanism of flower shape determination. Chrysanthemum morifolium, valued as an ornamental species given the diversity of its inflorescence form, is viewed as a model for understanding flower development in the Asteraceae. Yet, the underlying regulatory networks remain largely unexplored. Here, a transcriptomic survey of the Chrysanthemum morifolium variety 'Jinba' was undertaken to uncover the global gene expression profiles and identify the modules of co-transcribed genes associated with flower development. The weighted gene coexpression network analysis revealed important networks and hub genes including ray floret petals-specific coexpression network, disc floret petals-specific network, B and E class genes involved network and CYC2 genes network. Three ray floret petal-specific hub genes were also strongly transcribed in the ray florets of a selection of six diverse varieties and especially so in those which form ligulate ray floret petals. CmCYC2c was strongly transcribed in the distal and lateral regions of the ray floret petals, and also, along with CmCYC2d, in the tubular ray florets. Furthermore, CmOFP, belonging to the family of ovate proteins, was identified in the CYC2 genes network. CmOFP can interact with CmCYC2d that physically interact with CmCYC2c. This work provides important insights into the regulatory networks involved in flower development in chrysanthemum, and is informative with respect to the mechanistic basis of the regulation of flower shape.

Current achievements and future prospects in the genetic breeding of chrysanthemum: a review.

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a leading flower with applied value worldwide. Developing new chrysanthemum cultivars with novel characteristics such as new flower colors and shapes, plant architectures, flowering times, postharvest quality, and biotic and abiotic stress tolerance in a time- and cost-efficient manner is the ultimate goal for breeders. Various breeding strategies have been employed to improve the aforementioned traits, ranging from conventional techniques, including crossbreeding and mutation breeding, to a series of molecular breeding methods, including transgenic technology, genome editing, and marker-assisted selection (MAS). In addition, the recent extensive advances in high-throughput technologies, especially genomics, transcriptomics, proteomics, metabolomics, and microbiomics, which are collectively referred to as omics platforms, have led to the collection of substantial amounts of data. Integration of these omics data with phenotypic information will enable the identification of genes/pathways responsible for important traits. Several attempts have been made to use emerging molecular and omics methods with the aim of accelerating the breeding of chrysanthemum. However, applying the findings of such studies to practical chrysanthemum breeding remains a considerable challenge, primarily due to the high heterozygosity and polyploidy of the species. This review summarizes the recent achievements in conventional and modern molecular breeding methods and emerging omics technologies and discusses their future applications for improving the agronomic and horticultural characteristics of chrysanthemum.

Comprehensive characterization of a floral mutant reveals the mechanism of hooked petal morphogenesis in Chrysanthemum morifolium.

The diversity of form of the chrysanthemum flower makes this species an ideal model for studying petal morphogenesis, but as yet, the molecular mechanisms underlying petal shape development remain largely unexplored. Here, a floral mutant, which arose as a bud sport in a plant of the variety 'Anastasia Dark Green', and formed straight, rather than hooked petals, was subjected to both comparative morphological analysis and transcriptome profiling. The hooked petals only became discernible during a late stage of flower development. At the late stage of 'Anastasia Dark Green', genes related to chloroplast, hormone metabolism, cell wall and microtubules were active, as were cell division-promoting factors. Auxin concentration was significantly reduced, and a positive regulator of cell expansion was down-regulated. Two types of critical candidates, boundary genes and adaxial-abaxial regulators, were identified from 7937 differentially expressed genes in pairwise comparisons, which were up-regulated at the late stage in 'Anastasia Dark Green' and another two hooked varieties. Ectopic expression of a candidate abaxial gene, CmYAB1, in chrysanthemum led to changes in petal curvature and inflorescence morphology. Our findings provide new insights into the regulatory networks underlying chrysanthemum petal morphogenesis.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums.

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r 2 = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

Identification of favorable SNP alleles and candidate genes responsible for inflorescence-related traits via GWAS in chrysanthemum.

KEY MESSAGE: 81 SNPs were identified for three inflorescence-related traits, in which 15 were highly favorable. Two dCAPS markers were developed for future MAS breeding, and six candidate genes were predicted. Chrysanthemum is a leading ornamental species worldwide and demonstrates a wealth of morphological variation. Knowledge about the genetic basis of its phenotypic variation for key horticultural traits can contribute to its effective management and genetic improvement. In this study, we conducted a genome-wide association study (GWAS) based on two years of phenotype data and a set of 92,617 single nucleotide polymorphisms (SNPs) using a panel of 107 diverse cut chrysanthemums to dissect the genetic control of three inflorescence-related traits. A total of 81 SNPs were significantly associated with the three inflorescence-related traits (capitulum diameter, number of ray florets and flowering time) in at least one environment, with an individual allele explaining 22.72-38.67% of the phenotypic variation. Fifteen highly favorable alleles were identified for the three target traits by computing the phenotypic effect values for the stable associations detected in 2 year-long trials at each locus. Dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between highly favorable allele numbers and corresponding phenotypic performance were observed. Two highly favorable SNP alleles correlating to flowering time and capitulum diameter were converted to derived cleaved amplified polymorphic sequence (dCAPS) markers to facilitate future breeding. Finally, six putative candidate genes were identified that contribute to flowering time and capitulum diameter. These results serve as a foundation for analyzing the genetic mechanisms underlying important horticultural traits and provide valuable insights into molecular marker-assisted selection (MAS) in chrysanthemum breeding programs.

Overexpression of Phosphate Transporter Gene CmPht1;2 Facilitated Pi Uptake and Alternated the Metabolic Profiles of Chrysanthemum Under Phosphate Deficiency.

Low availability of phosphorus (P) in the soil is the principal limiting factor for the growth of cut chrysanthemum. Plant phosphate transporters (PTs) facilitate acquisition of inorganic phosphate (Pi) and its homeostasis within the plant. In the present study, CmPht1;2 of the Pht1 family was cloned from chrysanthemum. CmPht1;2 is composed of 12 transmembrane domains and localized to the plasma membrane. Expression of CmPht1;2 in roots was induced by Pi starvation. Chrysanthemum plants with overexpression of CmPht1;2 (Oe) showed higher Pi uptake, as compared to the wild type (WT), both under Pi-starvation and Pi-sufficient conditions, and also showed a higher root biomass compared to WT in the Pi-starvation conditions. Seven days after the P-deficiency treatment, 85 distinct analytes were identified in the roots and 27 in the shoots between the Oe1 plant and WT, in which sophorose, sorbitol (sugars), hydroxybutyric acid (organic acids), and ornithine (amino acid) of CmPht1;2 overexpressing chrysanthemum are specific responses to P-starvation.

Association analysis of drought tolerance in cut chrysanthemum (Chrysanthemum morifolium Ramat.) at seedling stage.

Understanding the genetic architecture is a prerequisite for crop improvement. The current research aimed to characterize the extent of genetic variation of drought tolerance harbored in a global collection of 159 chrysanthemum cultivars over 2 years. An average subordinate function value (ASFV), integrating the wilting index, the fresh weight retention rate, and the survival rate after re-watering under two drought-stressed trials, was used to quantify the level of drought tolerance. The performance of ASFV was generally correlated between the 2 years; and a high magnitude (0.95) of broad-sense heritability, coupled with the moderate genetic advance, was estimated for the ASFV. By applying MLM model with both population structure and kinship matrix as covariates association mapping identified 16 markers linked to drought tolerance, with the proportion of the phenotypic variation explained by an individual marker ranging from 4.4 to 7.6%. Of the eight markers predictive across the 2 years, four (E11M24-9, E3M2-8, E1M5-5, and EST-SSR34-3) were identified as favorable alleles for drought tolerance. Several cultivars that carry at least three of the four favorable alleles were identified as potential donor parents for future improvement of the drought tolerance. The findings provide an insight into the genetic basis of the drought tolerance in chrysanthemum and will, therefore, aid in developing new cultivars with enhanced tolerance against drought stress.

Dynamic and epistatic QTL mapping reveals the complex genetic architecture of waterlogging tolerance in chrysanthemum.

MAIN CONCLUSION: 37 unconditional QTLs, 51 conditional QTLs and considerable epistatic QTLs were detected for waterlogging tolerance, and six favourable combinations were selected accelerating the possible application of MAS in chrysanthemum breeding. Chrysanthemum is seriously impacted by soil waterlogging. To determine the genetic characteristics of waterlogging tolerance (WAT) in chrysanthemum, a population of 162 F1 lines was used to construct a genetic map to identify the dynamic and epistatic quantitative trait loci (QTLs) for four WAT traits: wilting index (WI), dead leaf ratio (DLR), chlorosis score (Score) and membership function value of waterlogging (MFVW). The h B2 for the WAT traits ranged from 0.49 to 0.64, and transgressive segregation was observed in both directions. A total of 37 unconditional consensus QTLs with 5.81-18.21% phenotypic variation explanation (PVE) and 51 conditional consensus QTLs with 5.90-24.56% PVE were detected. Interestingly, three unconditional consensus QTLs were consistently identified across different stages, whereas no conditional consensus QTLs were consistently expressed. In addition, considerable epistatic QTLs, all with PVE values ranging from 0.01 to 8.87%, were detected by a joint analysis of WAT phenotypes. These results illustrated that the QTLs (genes) controlling WAT were environmentally dependent and selectively expressed at different times and indicated that both additive and epistatic effects underlie the inheritance of WAT in chrysanthemum. The findings of the current study provide insights into the complex genetic architecture of WAT, and the identification of favourable alleles represents an important step towards the application of molecular marker-assisted selection (MAS) and QTL pyramiding in chrysanthemum WAT breeding programmes.

Genetic variation and association mapping of waterlogging tolerance in chrysanthemum.

MAIN CONCLUSION: Forty-five molecular markers were detected significantly associated with chrysanthemum' waterlogging tolerance, and four favorable parental lines were identified as potential donors for improving waterlogging tolerance in chrysanthemum. The productivity of chrysanthemum is downgraded by waterlogging soils, which has driven a search for germplasm showing an enhanced level of waterlogging tolerance (WT). As yet little is known regarding the mode of inheritance of WT in chrysanthemum. The study set out to characterize the extent of genetic variation for WT represented in a collection of one hundred chrysanthemum accessions by testing them under both greenhouse and field conditions. A membership function value of waterlogging (MFVW), which integrated a wilting index, a chlorosis score and the proportion of dead leaf in waterlogged plants, was used as a measure of WT. The variation for MFVW among plants grown in the greenhouse (two experiments) was generally higher than that generated in field-grown (one experiment) plants. The MFVW broad sense heritability was 0.82, and the phenotypic coefficient of variation (31.8 %) was larger than the genetic one (28.8 %). Association mapping (AM) identified 45 markers related to WT: 25 by applying the general linear model (GLM) + principal component (PC) model, 16 by applying the mixed linear model (MLM), 31 by applying the MLM + Q matrix model and 12 by applying the MLM + PC model. Of the associated markers, eight and two were predictive in two and three experiments within all models, respectively; the proportion of the phenotypic variance explained by the eight associations ranged from 6.3 to 16.4 %. On the basis of their harboring all four of the leading markers E2M16-2, SSR150-6, E19M16-1 and E10M10-12, the varieties 'Nannong Xuefeng', 'Qx097', 'Nannong Xunzhang' and 'Finch' were identified as potential donors for future improvement of WT in chrysanthemum.

Genetic diversity, population structure and association analysis in cut chrysanthemum (Chrysanthemum morifolium Ramat.).

Characterizing the genetic diversity present in a working set of plant germplasm can contribute to its effective management and genetic improvement. The cut flower chrysanthemum (Chrysanthemum morifolium Ramat.) is an economically important ornamental species. With the repeated germplasm exchange and intensive breeding activities, it remains a major task in genetic research. The purpose of the present study was to characterize the genetic diversity and the population structure of a worldwide collection of 159 varieties, and to apply an association mapping approach to identify DNA-based markers linked to five plant architecture traits and six inflorescence traits. The genotyping demonstrated that there was no lack of genetic diversity in the collection and that pair-wise kinship values were relatively low. The clustering based on a Bayesian model of population structure did not reflect known variation in either provenance or inflorescence type. A principal coordinate analysis was, however, able to discriminate most of the varieties according to both of these criteria. About 1 in 100 marker pairs exhibited a degree of linkage disequilibrium. The association analysis identified a number of markers putatively linked to one or more of the traits. Some of these associations were robust over two seasons. The findings provide an in-depth understanding of genetic diversity and population structure present in cut flower chrysanthemum varieties, and an insight into the genetic control of plant architecture and inflorescence-related traits.